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'T? Abstract 
Qh 

The cosmic-ray excess observed by PAMELA in the positron fraction and by FERMI 
and HESS in e~ + can be interpreted in terms of DM annihilations or decays 
into leptonic final states. Final states into r's or 4/i give the best fit to the ex- 
CN cess. However, in the annihilation scenario, they are incompatible with photon and 

neutrino constraints, unless DM has a quasi-constant density profile. Final states 
involving e's are less constrained but poorly fit the excess, unless hidden sector ra- 
O diation makes their energy spectrum smoother, allowing a fit to all the data with a 

combination of leptonic modes. In general, DM lighter than about a TeV cannot fit 
^ the excesses, so PAMELA should find a greater positron fraction at higher energies. 

O The DM interpretation can be tested by FERMI 7 observations above 10 GeV: if 

^ ^ the excess is everywhere in the DM halo, inverse Compton scattering on ambient 

^ light produces a well-predicted 7 excess that FERMI should soon detect. 
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1 Introduction 

Recently the PAMELA experiment observed an unexpected rise with energy of the j (e^ + e~) 
fraction in cosmic rays, suggesting the existence of a new component. The sharp rise suggests 
that the new component may be visible in the e~ + e+ spectrum. As worked out in [1] and 
further stressed in [2] , FERMI provides the first precise measurement of the e+ + e~ spectrum. 
The main purpose of this article is to analyze the implications of the first FERMI + e~ 
data [3], in conjunction with the new HESS measurements [4]. 

Although the peak hinted by previous data [5] is not confirmed, the FERMI and HESS [4] 
observations still demonstrate a deviation from the naive power-law spectrum, indicating an 
anomalous excess compared to conventional background predictions of cosmic ray fluxes at the 
Earth. Still, it is important to note that predictions of cosmic ray fluxes are model dependent 
and fraught with large uncertainties, therefore gauging the exact signiflcance or the correct 
interpretation is quite difficult. Astrophysical explanations have been put forth in an attempt 
to explain the data [6]: the most promising candidates are pulsars, which produce an excess 
with an energy spectrum that ranges from smooth to peaked to rastered. Thereby the FERMI 
measurement of the energy spectrum does not clarify the nature of the excess, and the most 
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interesting possibility remains open: the excess could be the first manifestation of Dark Matter, 
rather than a new background to Dark Matter searches. 

In this paper we concentrate on the Dark Matter (DM) interpretation. As we show below, 
the PAMELA, FERMI and HESS observations still point towards non-conventional DM models, 
with leptophilic final states, large "event rates" in our galaxy and a high mass scale ((9(TeV)). 
The reason for this is traced back to the fact that the new FERMI data does not show any 
sharp feature at low scale. Such a feature must be visible if the DM mass is low. Therefore, 
DM lighter than about a TeV is now excluded as an interpretation of the PAMELA excess. 
This conclusion further implies that the positron fraction should plateau or continue rising, but 
cannot go down at the higher energies being probed by PAMELA. Thus the DM scenario may 
be excluded in the future by PAMELA. 

Many models with these properties have been put forth to explain the excesses, and the range 
of possibilities span various annihilating and decaying dark matter scenarios. DM annihilations 
need to be enhanced with respect to what naively is suggested by thermal freeze-out in standard 
cosmology. This can be achieved either through the Sommerfeld effect [7, 2, 8], some resonance 
[2, 9], DM sub-clumps [10], or, of course, by a non-standard cosmology. It is therefore important 
to not only distinguish between the DM scenario and other astrophysical ones, but also to 
differentiate between the various DM models themselves. 

Indeed, while many models can explain the PAMELA anomaly, a large fraction of the 
existing models are inconsistent with other astrophysical measurements. Previous studies [2, 
11, 12, 13, 14] have analyzed some of the constraints and here we revisit the analysis once more 
in light of the new results. In particular, we consider the following set of measurements, 

• PAMELA positron fraction [15], anti-proton fraction [16], and preliminary electron flux [17]. 

• FERMI [3] and HESS [4] e+ + e" flux. 

• FERMI preliminary diffuse 7 rays [18]. 

• HESS and VERITAS observations of 7 rays from the Galactic Center [19] (including 
preliminary new HESS data), the Galactic Ridge [20] and dwarf galaxies [21, 22]. 

• SuperKamiokande up-going neutrino-induced /j,^ flux [23, 24]. 

The measurements above are insufficient to verify the DM picture. Ideally what one would 
hope for is some positive signal in a direct detection experiment. Nevertheless, there may still be 
important indirect evidence for DM by future measurements involving photons. These channels 
may shed light on whether the DM story is fully self consistent or not. In particular FERMI will 
measure additional diffuse gamma ray fluxes in our galaxy. This is an important measurement 
for two reasons. First, one expects DM to typically clump into subhalos within our Galaxy 
(as shown by A^-body simulations and theoretical arguments). These subhalos in turn can 
be a source of DM annihilation into channels with associated gamma rays and would provide 
an unambiguous signature for DM [10]. Second, barring the model dependency of including 
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DM with M = 3. TeV that annihilates into 4/i with crv = l.lx 10 cm^/s 




Positron energy in GeV Energy in GeV Photon energy in GeV 



DM with M = 3. TeV that annihilates into t^t with crv = 1.8 x 10 cm^/s 




Positron energy in GeV Energy in GeV Photon energy in GeV 



Figure 1: Sample DM fits. In the upper (lower) row we consider DM annihilations into 
fi~^ fi'' fi^ fi" (t^t^) with MED diffusion [25] and the NFW (isothermal) DM profile: all good 
fits are very similar. Left: the positron fraction compared with the PAMELA excess. Middle: 
the e+ + flux compared with the FERMI and HESS data. Right: the DM contribution to 
the diffuse photon energy spectra produced by bremsstrahlung (dashed red curve) and IC (black 
thick line); we also separately show the 3 IC components from star-light (red), CMB (green), 
dust (blue). 
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subhalos, the general DM profile extends above the galactic disk and should produce photon 
signals from Final State Radiation (FSR), Inverse Compton Scattering (ICS) and synchrotron. 
If DM exists in the form suggested by the charged particle fluxes, the signal from ICS in the 
diffuse gamma ray background as measured by FERMI should be an unambiguous clue as to 
whether DM is the cause of the excesses. 

In this paper we study DM that is consistent with all data including the preliminary 
PAMELA electron measurement, most recent FERMI electron plus positron measurements 
and the preliminary diffuse 7 measurements. We then make predictions for the additional 
FERMI diffuse 7 ray measurements and the future PAMELA positron fraction. In fig. 1 we 
demonstrate a sample good fit of the PAMELA and FERMI data and its prediction for the ICS 
spectrum that is an outcome of our analysis. The ICS flux in this example is representative of 
most models that fit PAMELA, FERMI and HESS. 

This paper is organized as follows. In Section 2 we describe the various final states in our 
DM calculations. We relate these final states to various models of DM including annihilating, 
decaying, and long lived intermediate particles. In Section 3 we outline the calculation for the 
contribution to the diffuse gamma ray background from ICS and synchrotron. We point out 
certain relations between the calculation of ICS and diffusion and bring up certain issues not 
discussed previously in the literature. We discuss a semi-analytic approximation to Inverse 
Compton scattering that illuminates some basic physical features. The reader might want to 
jump directly to Section 4, where we fit the various models to the existing cosmic ray data, 
identifying the parameter space viable for explaining PAMELA after the recent FERMI and 
HESS results. In view of the strong restrictions we make precise predictions for the upcoming 
diffuse gamma ray measurement by FERMI and comment on whether this will be a sufficient 
measurement to test the DM explanation. In the appendices, we describe some of the details 
of our approach to calculating neutrino fluxes and fitting procedure. 

2 Dark Matter models 

In this section we outline the "model" space that we cover in our study. Despite the existence 
of a cornucopia of models [8, 26, 27, 28, 29, 30, 31], we attempt to be as complete as possible 
in our scan while taking into account known constraints [2, 11, 12, 13]. By and large, DM 
models bifurcate into annihilating or decaying DM. In order to further explain the significant 
electronic activity and absence of hadronic activity, one either assumes a symmetry or otherwise 
the production of light states which are kinematically forbidden to decay into hadrons. In the 
latter case, the number of SM final states must be four or more. From the constraints point of 
view, the number of final states dictates the hardness of the injection spectrum leading to the 
following breakdown of the space of models: 

• Annihilating dark matter 

— Final states with 2 SM particles 
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— Final states with > 4 SM particles 
• Decaying dark matter 

— Final states with 2 SM particles 

— Final states with > 4 SM particles 

The bifurcation into annihilating versus decaying DM has consequences for indirect detec- 
tion because fluxes depend on the DM density p, as and p for annihilating and decaying 
models respectively. This implies that annihilating DM models are more constrained by the 
non-observation of a 7 excess from the Galactic Center, where the DM density is large (this 
bound can be weakened by assuming that DM annihilates into a long-lived particle [32]). De- 
caying DM models [26] have a lot more freedom. The excesses can be accommodated by 
choosing the DM decay rate, which unlike the thermal DM annihilation rate, is not linked to 
cosmology. Unfortunately, there is no definitive signal that allows one to differentiate between 
the decaying and annihilating DM scenarios. With this in mind our scan of model space will 
focus more on annihilating DM, but nevertheless include bounds and predictions for decaying 
DM. The relation between the rates of the two follows from a straightforward substitution: 
p^{av)/2]VP o pF/M, where r = 1/F is the DM lifetime. 

For annihilating DM, final states with 2 SM particles such as W^W~ or e^e~ are found in 
more conventional models of WIMPS [27, 28]. These models were already highly constrained 
prior to the release of the new FERMI data [2] and additionally require a non-thermal cosmology 
or large boost factor. They typically lead to either sharp features in the + e~ flux or large 
antiproton fluxes. The large antiproton fluxes are incompatible with the PAMELA data unless 
the scale is pushed to (9(10) TeV or maybe, with extremal astrophysical parameters, below 
~ 200 GeV [2]. The sharp features in the e"*" + e~flux that were a benefit in fitting the ATIC [5] 
data are now disfavored as we show in Section 4. 

The > 4 SM particle final state is a very interesting possibility for explaining the recent 
excesses. As has been noted, one needs a leptophilic final state and a large event rate. This 
can be obtained assuming that DM interacts with a lighter particle with mass m^, such that 
a long range attractive force generates a Sommerfeld enhancement to DM annihilations at 
low velocities. This lighter particle can no longer be the W boson [2], as DM annihilations 
into W^W~ does not provide a viable interpretation of the HESS and FERMI data. One 
possibility is to introduce an ad hoc new particle, for example by gauging — Lj. [2], allowing 
to explain the leptonic activity. A more elegant model of Sommerfeld-enhanced leptophilic DM 
annihilations is obtained [8] by assuming that the new particle only interacts with DM and 
is lighter than the proton, so that kinematics allow the new particles to decay only into the 
lighter leptons and pions. 

There are two leading effective operators of how a light hidden gauge group could couple to 
the SM and thus allow annihilations into visible SM particles. These are 

eF^M^/^v and c/i^/idm (1) 
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where F^J^ is the field strength of the dark gauge group, h is the SM Higgs, and /idm is a Higgs 
of the dark gauge group. The first operator describes the kinetic mixing of electric charge 
in our sector with the hidden kinetic term, and the second operator represents mass mixing 
through the Higgs. While both terms allow for DM annihilations into the light dark state, 
the subsequent mixing of the DM light state with SM particles determines whether it couples 
through charge (kinetic mixing with the photon) or mass (mass mixing with the higgs). In 
either case < 1 GeV is needed to explain the leptophilic final states, but the branching 
fractions are different, as they are proportional to the square of either the electric charge or the 
final state mass respectively. 

DM annihilations into two light hidden sector particles produce at least 4 body final states 
of SM particles. Indeed, as pointed out in [12], if the light hidden sector particle is charged 
under a hidden gauge group that is nonabelian (as preferred for phenomenological reasons [8]), 
final state radiation in the hidden sector can produce more hidden sector particles which sub- 
sequently decay into > 4 SM final states and can significantly effect final energy distributions. 
Alternatively, if the hidden sector has several higgs or multiple light particles there could be 
a cascade decay producing > 4 SM final states. Cascade decays and showering will not give 
identical spectra [12, 13], but since there is a great deal of model dependency in cascade decays 
we will only analyze the case of FSR. This can be computed in terms of the unknown gauge 
coupling constant in the DM sector, q;dm; and we consider several values for it to demonstrate 
the effects of multiple final light states that eventually decay to the SM. 

To demonstrate those models that are compatible with FERMI and PAMELA we choose to 
simulate the final states of 2yU, 2r, 4e, 4/i, An^, 4r, and additionally shower to create final states 
with > 4 SM particles. Other final states have been considered in the past [2, 12] and were 
shown to be excluded. In section 4 we concentrate on the bounds for each channel assuming 
these final states as fully exclusive, but we will also provide the constraints for models where 
the branching fractions are given assuming couplings proportional to charge. Other models 
which can generate > 4 SM final states but with different BFs then those explicitly shown in 
our paper [30] have bounds that can be inferred from a linear combination of our exclusive 
channels. 

If the light particle can decay into r^r^, it can also decay into pp. However due to kinemat- 
ical reasons, neglecting final state radiation, all protons have energy E > mpM/m^, which, for 
sufficiently high DM mass, M , is large enough to push the signal to energies not yet probed by 
PAMELA [12]. A hidden sector with light particles in the few-GeV mass range that decay into 
multi-r modes is presently motivated by the CDF multi-muon anomaly [33]. Since light parti- 
cles cannot be much heavier than 2mr (we fix = 4 GeV), the r spectra are not broad, and 
consequently the resulting e^,7,z/ spectra are just slightly broader than in the corresponding 
2r cases with DM mass rescaled by 1/2. However, r gives a large 7 fiux (due to the tt" decay 
chain), so that the r cases are only compatible with bounds on astrophysical 7 fiuxes [12] if 
DM decays (rather than annihilates) or if DM has a quasi-constant density profile (disfavored 
by A^-body simulations). 
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3 Testing the Excess with Diffuse Gamma Rays 

So far observations indicated excesses in cosmic ray e^. Unfortunately, the directionality of 
is lost by the galactic magnetic fields making almost isotropic in the sky. Observations of 

excesses therefore hardly allow one to test if it is present only locally or everywhere in the 
Milky Way. To have an unambiguous determination of the source direction we would need to 
observe a related excess in neutrinos or photons. 

Electrons generated by DM lose essentially all their energy via Inverse Compton scattering, 
e^7 — )■ on ambient light with average energy ~ eV. Such scatterings give rise to 

photons with larger energy Ey ~ E^{Ee/'mf,Y ~ 10 GeV, which is in the energy range being 
probed by FERMI. As discussed below, this DM ICS 7 flux is only marginally affected by 
astrophysical and DM distribution uncertainties. The reasons for this can be traced back to 
two observations: (i) Far away from the Galactic Center, the DM uncertainties are relatively 
mild, (ii) As we will see in Section 4, all DM models that fit the data predict roughly the 
same spectrum, as it is now mostly fixed by the new measurements (given the new FERMI 
and HESS results). Thereby the DM ICS spectrum is well predicted. As already illustrated in 
fig. 1 it is not much below the first FERMI diffuse 7-ray data, released for energies < 10 GeV 
in a specific angular region. Therefore, if the excess is due to DM, FERMI is expected to 
observe an associated 7 excess which is not sensitive to the specific DM model or DM density 
profile. Whether such an excess is seen or not, will decisively implicate on the DM (or any other 
mechanisms that produces in a spherical region away from the galactic plane) interpretation 
of the measured excesses. Alternative scenarios involve generated locally (e.g. by a powerful 
pulsar) or along the galactic plane (e.g. by supernovae). 

The ICS signal was computed in previous works in relation to the PAMELA measurement 
by running the GALPROP code [54] with some given astrophysical parameters [34]. Here 
we perform our own computation discussed in the rest of this section, where we justify our 
statement that the DM ICS signal has minor astrophysical uncertainties. We bring up some 
details not discussed previously in the literature on this subject. 

3.1 Energy Losses 

Energy losses of in the galaxy are due to two competing processes: 

1) Synchrotron radiation in the galactic magnetic fields B (with energy density ub = 
in natural units). 

2) IC scatterings on ambient photons, with energy density u^. 
The energy loss rate is given by. 



- 6(r, z, E) = —E^ 



(2) 



dt 3ml L i=i 

where cJr = Snr'^/S with re = aeui/irie- The light profiles u^{r^z) are relatively well known 
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Figure 2: Left: Energy spectra of the 3 galactic light components [55], normalized to unity. 
Right: The functions Ri{Ef.) which encode the relativistic corrections to the ICS energy loss. 
The black think line shows the function R{Ee) defined below Eq. (11). 



being composed of three components (hence the sum over i in the above equation): 

1. A constant CMB with energy spectrum du/dE.y = E^/n'^ /(e^''^'^ — 1) and mean energy 
(E^cmb) = 0.6 meV. 

2. Star-hght concentrated in the galactic center, with optical mean energy (E^star) = 1-1 eV. 

3. Star-light re-scattered by dust, with mean energy (-E^dust) = 0.01 eV. 

We plot the spectral shape of the three components in Fig. 2a. The functions Ri{Ee) in the 
above eq. (2) encode the relativistic correction to the non-relativistic Thompson limit of IC 
scattering and are plotted in Fig. 2b. They equal to unity at Ee <^ ml/ (Ei^) and exhibit a E'"^ 
suppression at higher energies. In practice, this relativistic effect must be taken into account 
for ICS on star-light for energies, Ee > 250 GeV. 

Magnetic fields, on the other hand, are not precisely known: they likely lie between 1 and 
10/iG, and may have the approximated profile [54], 

B{r,z) ^ ll/xG-exp(-r/10kpc- |z|/2kpc). (3) 

They give rise to energy losses into synchrotron radiation probed down to microwave frequencies. 
Thereby, while we know that essentially all the energy goes into photons, it is hard to 
establish the precise relative proportion of the two effects discussed above. Assuming the form 
of B{r,z) given in Eq. (3), Fig. 3 demonstrates the energy densities ub and u^i. We see 
that synchrotron energy losses are everywhere sub dominant by about one order of magnitude. 
In view of the relativistic effect discussed in the previous paragraph there is only one possible 
exception: synchrotron and ICS give comparable energy losses in the inner few kpc at E^ > TeV, 
where ICS energy losses on star-light are suppressed by the relativitic factor -Rstar- 
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Figure 3: Slices of the energy density profiles u{r, z) of star-light (red upper curve) re-scattered by 
dust (blue), CMB (green horizontal line), and of the presumed galactic magnetic fields (dashed). 
Left: Profiles as a function of the radius r, for a fixed z = 1 kpc. Right: Profiles as a function 
of the height, z, for a fixed r = 0. The values corresponding to B = 1 and 10/iG are indicated 
in the right plot. 



3.2 Computing 7 from Inverse Compton 

The IC process e^7 — )■ scattering of an with energy Eg and isotropic initial direction 
on isotropic 7 with energy gives 7' with the following 7' energy spectrum [35]: 

- 2vrr2-f-I/ic(g,e), (4) 



dE'^ dt ' El E^ 

/ic(g,e) = 2glng+(l + 2g)(l-g) + :^^^(l-g). (5) 

Here Ee ^ rrie and 6, F, g are the dimensionless variables defined as: 

E'^ AE,E, e 

Ee ml r(l — e) 

E^i lies in the range E^/ E^. < e < r/(l + r). The non-relativistic (Thompson) limit corresponds 
to r ^ 1, so that e ^ 1, the last term in /ic is negligible, and < g < 1. The total energy loss 
rate in eq. (2) is recovered by integrating eq. (4) over E'^. 

The ICS flux $^ = dN^/dSdt obtained for a given line of sight and generic energy and 
spatial distributions of initial and 7, is found to be, 

dEydQ ^ 2 ""'"Acs. ^JJ dEedK,^^^^"^' ^ ' 

Here we do not include the anisotropic correction to 7 due to the fact that light is emitted 
preferentially from the galactic plane. The effect in the FERMI region is (9(10%) [56]. The 
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DM ICS flux depends on the density / = due/dE. It is computed by numerically solving 
the diffusion-loss equation 

.K(E).vy-§-^Hrr^.m^\{ii)\^^)'-^ (8) 

where dNe/dE is the spectrum produced by one DM annihilation while b{r,z,E) is given in 
Eq. (2). For more details on the above diffusion-loss equation and its approximate solutions see 
e.g. [36, 37, 12] and references therein. On the practical level, we work on a grid in cylindrical 
coordinates r, z, E and discretize the differential equation (8) and its boundary conditions into a 
system of thousands of linear equations, easily solved with Mathematica. We take into account 
the spatial dependence of the energy loss term b assuming that each of the light components 
u^i, ub, can be factorized in space and energy as discussed below. 



3.3 Approximating the Diffuse 7-ray Flux 

In order to gain more intuition, it is beneficial to study the ICS spectrum in an a approxi- 
mate manner. Let us, therefore, obtain a simplified expression for the diffuse gamma-ray flux 
which illuminates the physics. The photon spectrum is the sum of 3 components (star-light, 
dust, CMB) which have position-dependent intensity, and roughly the same energy spectra at 
any point: dn^/dE^ = J2i f'yi{E^)u^i{r) / {E^i) where u-yi is the total energy density, and the 
functions f^i are plotted in fig. 2a and normalized to unity: / dE^f^i = 1. 

Assuming that diffusion is negligible^ the energy spectrum of generated by DM has 
the same shape at any position: 



dE dE AarUtot 2 \ M J ' ' Je dE^ 

where / is defined as in [36, 37]. Therefore the ICS spectrum simplifies to 

where the dimensionless factors Jic and Gic respectively encode astrophysics and particle- 
physics: 

We explicitly see that the overall amount of energy Utot = ub + J2iU-yi does not matter, but 
only the ratios among the different components, up to the relativistic correction R{Ee) = [ub + 
J2i UyiRi{Ee)]/utot- Unfortunately, in general this last factor is position-dependent and therefore 
Giic is not purely determined by particle-physics. Nevertheless, this latter dependence is weak 
and to leading order Guc encodes the particle physics information which enters the diffuse 



^In the halo-function formalism [.36, 37] this corresponds to assmning unity halo-fuction, / ~ 1. 
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region 




DM annihilation 




DM decay 




10° < \b\ < 20° 


NFW 


Einasto 


isothermal 


NFW 


Einasto isothermal 


<^star,IC 


4.6 


6.0 


2.9 


2.4 


2.6 


2.2 


<^dust,IC 


1.1 


1.4 


0.78 


0.67 


0.70 


0.63 


>/cMB,IC 


1.0 


1.2 


0.86 


1.4 


1.4 


1.5 


J 


7.4 


9.2 


5.0 


4.9 


5.0 


4.7 



Table 1: Astrophysical factors Jj for Inverse Compton within the 10° ^20° ' region observed 
by FERMI. Since the region is large and away from the GC, the Ji only have a minor dependence 
on the DM profile. 



spectrum. In fact, assuming that IC scattering on star-light dominates energy losses up to 
an energy E^, it is given by R{Ee) ~ i?star(niin(ii^e, -E*)), where astrophysics is all condensed 
in the parameter i?*, which presumably is large enough not to be crucial. The function R is 
plotted in fig. 2b for the FERMI region. 

It is interesting to note that in the case where synchrotron energy losses are subdominant 
with respect to the ICS energy losses, Jnc = J, where J is the usual factor that encodes 
the astrophysics of DM annihilations into photons (with energy spectrum dN^/dE^ per DM 
annihilation): 

d^. _/j^^dN^(Pg\\ j=[dn[ ±(pM]\ (12) 



dE^ Stt dE^ \MJ 

In table 1 we list the astrophysical J factors for the '10° 20°' region observed by FERMI. It 
is indeed apparent that J amounts to the bulk part of the contribution. Moreover, since the 
region probed is far from the GC, the various Ji factors only weakly depend on the DM profile 
and the predicted spectrum is therefore robust. 

As a word of caution we stress that the J integral should be computed only inside the 
diffusion volume. For the above table we assumed the MED propagation parameters {L = 
4kpc), and, in the case of the '10° ^ 20°' region this cut reduces J by about 10% with respect 
to its full volume value (L = oo). A significant reduction is found only for the implausible 
MIN configuration {L = Ikpc).^ Finally, throughout our analysis below, we used both the 
approximation discussed above and the exact solution to the diffussion-loss equation. We 
found that the approximation for the ICS fiuxes differs by a factor of 0{2). This is illustrated 
in Fig. 4. 

The plots in Fig. Ic show the result of our full ICS computation. We have separately shown 
the DM predictions for the three ICS components within the '10° -v- 20°' region so far observed 
by FERMI. The DM model plotted gives a best fit to the FERMI, PAMELA and HESS 
excess. In fact, as mentioned above, the electron spectra turn out to be rather similar for the 



^ We compute ICS neglecting the contribution of generated outside the diffusion cyUnder. Indeed, despite 
energy losses due to ICS, their mean free path is « 100 kpc for E w 1 TeV and u « eV/ cm'^, larger enough 
than the size of our galaxy, that their contributions to the ICS 7 flux gets mostly lost. 
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Figure 4: Inverse Compton: exact vs approximated. We compare our full calculation for 
ICS (solid curves) with the diffusion-less 1 = 1 approximation of eq. (10) (dashed curves) for 
the '10° 20° ' region and DM annihilating into with av = lO^^^cm^/s and NFW with 

MED propagation. Additionally we plot the individual contributions to ICS from CMB (blue), 
dust (green) and starlight (red), while the total contribution is in black. The ICS approximation 
is good within a factor of two. 

majority of models that fit FERMI, PAMELA and HESS and therefore the ICS spectra are 
very similar in the various cases. We provide more evidence to this statement in section 4. 
This observation is in contrast to the FSR 7 fluxes (red dashed lines in Fig. 1) which are highly 
model dependent [11, 12, 13] and typically dominate over ICS at higher energies close to the 
DM mass. 

In conclusion, the prediction for the ICS 7 flux is robust: in the '10 20' region the DM 
contribution should be visible at > 100 GeV, and other regions offer better sensitivities. 
Reducing the IC DM signal would require unexpected astrophysics: either very large galac- 
tic magnetic fields (such that the detectable signal anyhow moves from ICS to synchrotron 
radiation) or perhaps a very thin diffusion cylinder. 

4 Interpreting the PAMELA and FERMI Observations 

We now demonstrate the viable models that can fit all the available data for charged cosmic 
rays including the new FERMI and HESS data. For our analysis we fit to the PAMELA 
e"*"/ {e~^+e~) data and the e~^ + e~ data from FERMI and HESS. We also include the preliminary 
un-normalized PAMELA e~ spectrum available at [17], although it has a minor impact. We 
then describe the predictions and bounds for ICS as discussed in Section 3, and include bounds 
from several other observations previously studied in [11, 38, 37, 12]. These include the bounds 
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DM annihilation 



DM decay 




Figure 5: Global fit to PAMELA, FERMI and HESS data. The labels on each curve indicate 
the primary DM annihilation (left) or decay (right) channel. In the left panel a hypothetical flux 
from a pulsar is also plotted, with an assumption that the flux is given by ^ = E^^e~^^^'^ . In 
the left panel all final states for DM annihilation do not include hidden sector FSR, except for 
the curve labelled Afi — sh. This curve demonstrates that by including the hidden sector shower 
the is significantly improved and is as good of fit as any other hypothesis. 



from the HESS photon measurements in the Galactic Center and Galactic Ridge and from 
up-going muons measured at SuperKamiokande. 

The FERMI data is conservatively fitted adding in quadrature statistical with systematic 
uncertainties independently for each data-point. We consider uncertainties on the smooth DM 
halo profile p(r), propagation, and the spectral index and normalization of the e"*" and of 
the e~ astrophysical backgrounds. We keep fixed the local DM density, = 0.3 GeV/ cm^. 
Changing it would be equivalent to an overall rescaling of the DM annihilation or decay rate, 
which renormalizes in the same way all indirect DM observables. Therefore, the comparison 
between the regions favored by the excesses and the constraints from 7 and u observations 
remains fully meaningful. In Appendix B we describe further how our fit is performed, and 
here we simply summarize the main points. 

Other, less established bounds, are not included. In particular, ref. [39] finds that angular 
regions distinct from the one observed by FERMI and so far observed only by EGRET, provide 
stronger constraints. We do not use here the controversial EGRET observations. Once FERMI 
will present data corresponding to other regions, it will be easy to establish bounds with the use 
of our approximation described in the previous section. This is done by simply rescaling our 
predictions for the '10° -r- 20°' region using the new J factors for the additional (yet unknown) 
regions. Furthermore, we do not consider the 'WMAP haze' [40] which is a hint that a possible 
excess in synchrotron radiation could be due to DM. The haze has been shown to be consistent 
with a wide variety of DM masses and final states and therefore will not constrain the space 
of models compared to other measurements. It would however, be interesting to study in more 
detail the precise predictions for the haze for those models that can fit the rest of the data. 
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Fig. 5 shows the as function of the DM mass for various DM annihilation (left) or DM 
decay (right) modes. We find that, independently of the non-observation of an excess in the 
p PAMELA data, only some leptonic modes can reproduce all data. Here HESS observations 
play a key role demanding that the excess terminates in a sharper way than what typical 
of non-leptonic channels, irrespectively of the DM density profile. DM heavier than 10 TeV 
that annihilates or decays into light quarks still provides a reasonable fit to PAMELA and 
FERMI data, if fitted conservatively. However it is disfavored by the HESS e"*" + e~ data, and 
presumably the photon data as well ( in the annihilating case). 

The only spectral feature that can allow one to discriminate the various modes lies at the 
high end of the spectrum between 1 and 3 TeV, where we only have the electron HESS data 
which is less precise than the FERMI data. The leptonic channels can be ordered according to 
the sharpness of their end-point: 

2e > 4e ~ 2/i > 4^ ~ 2r > 4r. (13) 

The main difference with respect to the previous ATIC + data [5] is that a peak is no 
longer present, and therefore the 2e mode (namely, DM annihilations or decays into e~^e~) is now 
excluded independently of the photon bounds, since it predicts a too sharp end-point. All other 
leptonic modes we consider provide comparably good fits to the data, that therefore cannot 
discriminate which (combination of) modes is the correct one. We also explored the spectra 
produced by polarized fi or r, similarly finding that present data do not allow to discriminate 
the various possibilities. 

The FERMI data alone, becomes more constraining if one takes into account the appro- 
priate correlations among the systematic uncertainties at different data points. Assuming that 
systematics only affects the overall energy scale, and that our description of the background 
is still valid at the level of precision of FERMI statistical errors, data would imply a prefer- 
ence for the smoother 4/x or 2r modes. As stated above, assuming instead that FERMI data 
must be fitted conservatively (so that they are consistent with a power law with no spectral 
features) and dropping the HESS data (which implies a termination of the excess at ~ 2 TeV), 
non-leptonic modes are consistent with the data if DM is heavier than ~ 10 TeV and provide 
fits as good as the leptonic modes for M ~ 30 TeV. In particular. Minimal Dark Matter (which 
predicted the PAMELA excess as Sommerfeld-enhanced DM annihilations into W^W~' with 
M 9.6 TeV and no ATIC peak [37]) is no longer a viable DM interpretation of the + e~ 
excess, if HESS or FERMI conclusively establish that the excess terminates around 2 TeV (or 
if the DM profile is found to be sufficiently dense at the GC to be constrained by HESS). 

Therefore a future more precise measurement of the e"*" + e~ spectrum around the end-point 
of the excess at ~ 2 TeV will be very important to disentangle the DM annihilation modes. 
This is illustrated by the best-fits shown in fig. 8b for the 4£ annihilation modes. We here fixed 
MED propagation and the Einasto DM density profile, so that each best-fit could be further 
slightly improved by adjusting the propagation and the DM profile. We see that the various 
channels mildly differ at the (2-^3) TeV energy at which the + e~ excess terminates. We also 
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see that in all cases the other feature present in the FERMI data, namely the mild hardening of 
the e"*" + e~ spectral slope around 100 GeV, can be simply attributed to the takeover of the DM 
component, without indicating any special spectral feature in it. Finally, all best fits predict 
that the the positron fraction must continue to go up to about 1 TeV in an essentially unique 
way dictated by the FERMI data. 

4.1 Best e± Fits vs 7 and u Bounds 

Figures 6, 7 and 10 show the best-fit regions in the (M, (crv)) plane for DM annihilations 
and in the (M, r) plane for DM decays. Each panel assumes a specific mode and a specific 
DM profile. In each panel the red regions are favored by the global fit of FERMI, HESS and 
PAMELA data at 3 and 5cr (2 dof). The green bands are favored by PAMELA (at Sa for 1 
effective dof; we here used the lowest-energy FERMI bins to fix the overall flux, finding a {av) 
which is a few times lower than the best-fit bands of [11]. This could only be due to distinct 
theoretical computations of the fluxes and the fit to the ATIC peak). 

The plots show that the best-fit DM mass and annihilation cross section or decay rate is very 
different in the various cases, ranging from 1 TeV (for the fi~^fj^~ channel) to almost 10 TeV (for 
the 4r channel). Thereby the channel could be discriminated looking at different observables. 
In view of the multi-TeV DM masses, it will be difficult to test these DM scenarios at LHC 

The most promising observables are the 7 and u fluxes generated by DM annihilations or 
decays. As all data for the moment shows no excess, we compare the regions that best fit the 

excesses with the regions excluded by other observations. Even if the best-fit regions 
are almost the same irrespectively of the DM density profile, such comparisons must be done 
for a given DM density profiles. Indeed, photon measurements point to their source, and so 
the associated bounds strongly depend on the DM density profile. We consider the NEW [41], 
Einasto [42] and isothermal [43] DM density profiles. 

The shaded regions in the figures are excluded by the various 7 and u constraints: 

1. The blue continuous curves labeled as 'GC-7' shows the bound obtained by imposing that 
the HESS observations of 7 rays from the Galactic Center region (with 0.1° angular 
size) are not exceeded at 3a in any data point. We updated the Galactic Center constraint 
from [11] in light of the recent preliminary data from HESS [19], shown in fig. 8d. 

2. The blue dot-dashed curves from [11] labeled as 'GR-7' show the corresponding constraint 
from 7 observations of the Galactic Ridge region [20] . In view of its larger angular 
size (a rectangle of size 0.6° x 1.6°) these bounds have a milder dependence on the unknown 
DM profile. Again we conservatively imposed that the DM contribution does not exceed 
any data point at more than 3cr. 

3. The blue dashed curves labeled as 'dS-7' show the best bound from non-observation of 
7 from dwarf Spheroidals. We proceed as in [11] where the HESS observations of 
these objects were first studied to test the DM interpretations of the PAMELA excess. 
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Figure 6: Direct DM annihilation. We compare the region favored by PAMELA (green 
bands) and by PAMELA, FERMI and HESS observations (red ellipeses) with HESS observations 
of the Galatic Center [19] (blue continuous line), Galactic Ridge [20] (blue dot-dashed), and 
spherical dwarfes [21, 22] (blue dashed), FERMI observations in the 10° -i- 20° ' region and of 
observations of the Galactic Center at radio-frequencies v = 408 GHz [44] (dashed red lines) 
and atv lO"*^^ Hz by VLT [45] (upper purple lines, when present, for equipartition and constant 
magnetic field) . See discussion in the text for remarks regarding the validity of the constraints. 
We considered DM annihilations into e^e~ (left column), fi^fj^" (middle), t^t^ (right), unity 
boost and Sommerfeld factors and the NEW (u^per row), Einasto (middle), isothermal (lower) 
DM density profiles in the Milky Way. 
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Figure 7: One-step DM annihilation. ^45 in fig. 6, here for DM annihilations into 4e (left 
column), 4/i (middle), 4r (right), via a light intermediate new particle. 
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Spheroidal DM annihilation DM decay bound on the 7 flux 

dwarf logio(i2 in GeV^cm"^) logio(vC in GeVcm"^) in cm~^s~^ 

Sagittarius 19.4 ± 1.0 18.7 ±0.9 <^.,{E^ > 250 GeV) < 3.6 x lO^^^ [21] 

Draco 18.6 ± 0.4 17.5 ± 0.1 <^j{E^ > 200 GeV) < 2.4 x 10"^^ [22] 

Ursa Minor 18.8 ± 0.8 17.6 ± 0.2 (^^{E^ > 200 GeV) < 2.4 x 10"^^ [22] 

Willman 1 19.6 ± 0.6 17.5 ± 0.5 <^j{E^ > 200 GeV) < 2.4 x lO^^^ [22] 

Table 2: Central values and Icr errors for the astrophysical factors C defined in eq. (14) for 
DM annihilations and decays in Milky Way dwarf galaxies as compiled in [46]. Note that the 
bounds coming from Sagittarius should be taken with extreme caution since it is being visibly 
tidally disrupted by our galaxy, and therefore the DM profile or even ascribing that the dwarf 
has any DM is potentially suspect. 



but now adopting the DM luminosities C and the corresponding errors estimated in [46] 
and defined such that 

d^^ C (av) dN^ , , 

— - = — ^ — -. — - (14) 
dE^ 47r2M2 dE^ ^ ^ 

For DM annihilations, in the 'NFW panels we assume the central values of £; in the 
'isoT' panels we assume the 3a lower limits on C and in the 'Einasto' panels we assume 
the 90% CL. lower limits on C as in [46]. For DM decays, we will in all cases consider 
the 3a lower limits. We stress that the bounds coming from the Sagittarius Spheroidal 
dwarf should be taken with caution since its close proximity to the Milky Way causes it 
to be tidally disrupted which renders the use of conventional DM profiles questionable. 
This is partially refiected by the larger uncertainties in table 2, such that Sagittarius does 
not provide the dominant constraint. 

4. Radio observation of the inner 4 arc seconds of the Galactic Center constrain syn- 
chrotron radiation produced by in the local magnetic fields B. While the uncertainty 
on B seems not important (as anyhow loose their energy via synchrotron radiation), 
the local DM density is very uncertain, so that the resulting bound on DM annihilations 
(red dashed curve labeled as 'GC-radio', from [11]) is very strong (negligibly weak) if a 
NFW (isothermal) profile is assumed. We plot this bound in light colors because it relies 
on the untested extrapolation that the growth of the DM density profile as r — )■ holds 
down to r ~ pc. 

5. The black continuous curves labeled as V show the bounds on neutrinos from the 
Galactic Center regions, as observed by SuperKamiokande in regions with angular 
size between 3° and 30°: in view of the large observed region such bounds have a weak 
dependence on the uncertain DM density profile. We consider the through-going muon 
fiux and impose that the contribution due to DM does not exceed the observed rate at 
3a. In appendix A we describe how we compute this effect. As neutrino detection cross 
sections increase with energy, these constraints are significant at larger DM masses. 
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6. Finally, the red continuous curves labeled as 'IC are relative to the preliminary FERMI 
observation in the '10°-i-20°' region of 7. As discussed in section 3, in view of the low 
energy measurement, the DM signal is dominated by ICS 7 rather than by Final State 
Radiation. 

Below we summarize the general lessons learned. 
4.2 DM Annihilations 

Concerning DM annihilations, 2e, 2/i, 2r, 4r modes are strongly constrained and are compatible 
with 7 bounds only if DM has an effective density profile that does not significantly grow 
towards the Galactic Center: either the isothermal profile, or any profile in models where DM 
annihilations proceed via extra long-lived particles [32]. 

Only the 4e and 4/i modes are (marginally) compatible with the Galactic Ridge HESS 
observations also for a NFW or Einasto profile. Even these cases can already be strongly 
disfavored by performing a less conservative global fit of all HESS data points in terms of an 
astrophysical background plus the DM excess. This conclusion is changed if hidden sector 
showering plays a role, as we discuss below. Going from 2£ to M (where £ = e or /x) the FSR 7 
flux is changed in shape [12] and normalization [12, 11]. However, it is difficult to further lower 
the 7 constraint below what shown in the 4/i panel. 

The expected neutrino signal is comparable to the present bound on neutrinos from the GC, 
again favoring a quasi-flat DM density profile or the 4e mode that does not give any neutrinos. 
Present SK bounds are limited by statistics but also by the z/ atmospheric background, so that 
an improvement would need selecting a sample of higher energy neutrino-induced through- 
going muons. This was already attempted in SK, that was able to select a sub-sample of 
muon- showering events, namely of muons with > TeV energy in the detector such that muon 
energy losses are dominated by radiative processes [24] (as here discussed in eq. (15)). Such 
sub-sample is expected to imply bounds comparable to the bounds considered here. 

4.2.1 Hidden Sector Shower 

As discussed in Section 2, an interesting possibility for models that contain a light vector that 
the DM can annihilate through is a hidden sector shower. Showering has two primary effects: 
(i) it softens the spectrum and (ii) it increases the multiplicity. As we see in Figures 5 and 9, 
the hidden sector shower has a nontrivial effect, and improves the of our fits. The hidden 
sector shower helps reduce the by reducing the steepness of the + e~flux and better fitting 
the FERMI electrons. However, as we can see from Fig. 9, showering does not improve the 
for all values of odm- For smaller values of odm the spectrum of visible particles is softened 
without increasing the multiplicity of final state particles greatly. On the other hand, for larger 
values of odm, the bounds will strengthen as a function of odm as the multiplicity becomes too 
large, generating too many 7's. 
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Figure 8: Best fits for DM annihilations. We assume the MED propagation model and the 
Einasto profile with a = 0.17. We plot the best fit models for DM annihilations into Ac, 4/i, 4r, 
and final states dictated by coupling through kinetic mixing with the photon (m^ = 650 MeV) 
see figure 9. All curves include showering with udm ~ 0.1 which increases the goodness of their 
fit compared to the unshowered spectra. The plots in the upper row from left to right are for the 
PAMELA positron fraction, + e~ flux recently measured by FERMI and HESS, and the ICS 
+ FSR predictions for these models. In the bottom row we plot from left to right, the photon 
predictions for the HESS measurement of the Galactic Center, Calactic Ridge, and the bounds 
coming from SuperK for those models which create v 's. 



4.2.2 A Quasi-constant DM Density and Long Lived Intermediate States 

DM annihilations into 2/i and 2r are still compatible with bounds on the associated 7 flux if the 
DM density does not significantly grow towards the Galactic Center. This possibility is realized 
in practice by plotting the 'isothermal core' DM density profile. This profile is disfavored by 
A^-body simulations and it has no a priori theoretical motivation. However, in principle there 
could be some weakening of the more cuspy DM profiles preferred by A^-body simulations 
when baryons, which should be important near the center of the galaxy, are included in future 
simulations. To explore the effects of a shallower DM density profile we have shown in the 
lower row of Fig. 1 the best fits for DM annihilations into 2r assuming the extremal isothermal 
profile and MED propagation. 

For the case of > 4 SM final states there is one other possibility to obtain an effective quasi- 



21 



2Q I I , I , I , , I I ^ JO ^ — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ^ 

1000 1500 2000 3000 5000 7000 10000 0.00 0.02 0.04 0.06 0.08 O.IO 

DM mass in GeV 0!dm 

Figure 9: Hidden sector shower. Left: the with and without shower, for various light 
gauge boson masses of 250, 450, and 650 MeV (blue, red, green respectively) that kinetically 
mix with the photon. Additionally, a curve for a 4r final state without shower is shown for 
comparison of the overall x^- Right: Contour lines for various bounds in the a vs a^M plane, 
for a DM mass of 3 TeV and light gauge boson of mass 650 MeV with an Einasto profile. 

constant DM density profiles. DM can annihilate through a light state that can be long lived on 
astrophysical scales [32]. There is one additional parameter, the lifetime c?, that effectively acts 
such that indirect DM signals must be computed replacing the DM density profile p(r) with its 
value averaged over a length d. \i d ^ Tq a, constant p is effectively obtained. In our analysis 
this possibility is effectively described by the quasi-constant isothermal profile. For smaller 
values of d the effective p(r) flattens to a constant only at r < d rather than diverging as in 
the NFW or Einasto profiles; the main consequence is relaxing the bounds from 7 observations 
from the GC region at r < li that we consider. 

4.3 DM Decays 

The PAMELA and FERMI excesses can be explained if DM decays into leptons with a life-time 
about 10^ times longer than the age of the universe, which is the typical lifetime of a TeV-scale 
particle that decays via a dimension 6 operator suppressed by the GUT scale, 10^^ GeV [47, 48]. 
As discussed in [47] this scenario is compatible with the 7 bounds even for a NFW scenario, 
as the DM decay rate is proportional to the DM density p, while the annihilation rate is 
proportional to and thereby significantly enhanced close to the Galactic Center. 

Here we remark that cosmology offers other constraints on DM annihilations: BBN [49], 
extra-galactic 7 [52], gas heating [50] and reionization [50, 51]. These constraints are significant 
thanks to the fact that the DM annihilation rate was enhanced in the early universe by the 
square of the larger DM density. The latter constraint is solid, while the first three ones are 
uncertain, as they depend on the unprecisely known structure formation history at redshift 
z < 100. These constraints are less significant in the alternative interpretation in terms of DM 
decays. 

We here assume that DM decays with 1/2 branching ratio in leptons and in anti-leptons. 
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Figure 10: DM decay. As in fig. 6, here for DM decaying into ^'^ (middle), t^t^ (right), 
4/i (left). We do not consider decay modes into e^e~ , as they do not allow to fit the FERMI 
data. 
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Figure 11: Best fits for DM decays. As in fig. 1, now for DM decays into 4/i (red curve), 
2t (blue curve) assuming the MED propagation model and the NFW profile. 



although DM present thanks to a baryon-hke asymmetry could decay only into anti-leptons. 

Fig. 5b shows that DM decays can fit the PAMELA, FERMI and HESS excesses as 
well as DM decays, with some minor differences, mostly due to the fact that only in the DM 
annihilation case a sizable amount of lower-energy can reach us from the Galactic Center, 
giving rise to a smoother energy spectrum. Decay modes of fermionic DM into W^fi^ or 
W^T^ provide good fits to observations, but together with a p excess which is strongly 
disfavored by PAMELA p observations [2]. 

Fig. 10 and 11 show that DM interpretations of the PAMELA, FERMI and HESS 
data in terms of leptonic DM decays are compatible with all the constraints we considered, 
the strongest one being relative to ICS from the '10° ^ 20°' region observed by FERMI. Other 
regions (so far observed only by EGRET, with a problematic energy calibration) offer possibly 
stronger constraints. In the near future FERMI will release data in other regions that can be 
used to constrain this scenario. As discussed in section 3, one can easily extrapolate the fluxes 
to other regions using the approximations we have employed. 

4.4 Light Dark Matter and the Electron Component 

DM lighter than about a TeV is firmly excluded in a model-independent way as an interpretation 
of the PAMELA excess. Indeed according to PAMELA, the positron fraction reaches ~ 15%, 
so that a drop of at least > 20% should have been present in the e+ + e~ spectrum at an energy 
just below the DM mass. This is not seen in the FERMI data, which have a ~ 1% statistical 
uncertainty together with a ~ 5% systematic uncertainty that cannot fake this drop. This 
conclusion was already derived in [2] in the light of the less precise ATIC data, confirmed at 
lower energies by FERMI. 

Similarly, the lack of an edge in the + e~ spectrum implies that DM that dominantly 
annihilates or decays into 2e is now firmly excluded, and one can constrain the BR of the 
subdominant primary channel. Assuming that DM annihilates into leptons with lepton 
flavor components BRe, BR^, BR^ that sum to unity, BRg + BR^ + BR,- = 1 one can constrain 
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Figure 12: Fits for the lepton flavor component in DM annihilations. The contour lines 
correspond to 90%, 99%, 99.9% confidence levels for 2 d.o.f. In the left panel we consider an- 
nihilation into 4i modes with Einasto profile. In the center we show 2i modes with isothermal 
profile, and on the right we show decays into 2i modes with isothermal profile. Pure flavors are 
obtained at the corners of the triangles, and the distances BR^ from the opposite sides (which 
sum up to unity) indicate the flavor proportion at each point. The mostly-r region below the 
dashed (solid thick) curves is excluded by the Galactic Ridge (v) bounds. 

the allowed regions in light of the PAMELA, FERMI and HESS data together with the other 
constraints. To clarify the situation we show the constraints on the branching fractions for 
three distinct cases: 

1. DM annihilating into Ai with Einasto profile. HESS and SK provide strong constraints 
on the associated 7 and u fluxes, which must be taken into account. 

2. DM annihilating into 2i with isothermal DM profile. Steeper profiles are excluded by 
HESS and we therefore restrict ourselves to this profile only. 

3. DM decaying into 2i with NEW profile. As shown above, decaying DM evades all bounds 
and is therefore viable for all DM density profiles. 

Fig. 12 shows the fit for these cases, marginalized over the DM mass and other parameters. 
The FERMI data has been fitted conservatively. Nonetheless we find that FERMI disfavors 
DM that annihilates democratically into all with branching ratios equal to 1/3 (point in 
the center of the triangle). This is the situation in many flavor-universal models. A similar 
conclusion holds for DM decays into W^i^. On the contrary the 4e mode is more compatible 
with data, as it gives a smoother termination of the excess at ~ M. 

5 Conclusions 

We explored how DM can interpret the excesses in cosmic rays measured by FERMI, 
PAMELA and HESS. The results of the HESS and FERMI data indicate two specific things 
for any DM interpretation. First, the spectrum becomes steeper at around 1 TeV implying 
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the scale of DM needs to be around a TeV. Second, the FERMI data is smooth from 100 GeV 
to around a TeV implying that a light DM explanation of PAMELA is inconsistent with the 
FERMI data. 

We have found that the PAMELA, FERMI and HESS data can be interpreted in terms of 
admixtures of leptonic final states for DM annihilations or decays. We considered the e~^e~ , 
yU^/i^, r+r^ 4e, 4/i, Air and 4r modes. 

• The e"*"e~ mode is excluded. 

• Modes involving r provide one of the best fits to the spectra. However the r decays into 
7r° imply a large 7 flux, so that r modes are compatible with bounds from 7 observations 
only if a) DM has a quasi- const ant density proflle, disfavored by A^-body simulations; or 
b) DM annihilates into a particle with an astro-physically long life-time; or c) DM decays. 

• Annihilations into 4yU provide a good fit to the spectra (for M ^ 3 TeV) and are 
marginally compatible with 7 and u bounds for the Einasto profile favored by A^-body 
simulations. Again, decays are cleanly compatible with all constraints. 

• The 4e mode provides a poorer fit to the FERMI spectrum. 

• Annihilations into 4 body final states that are dictated by a hidden sector coupling to the 
SM proportional to charge are marginally consistent with all data. 

• Models that have a hidden sector shower can provide a better fit, by smoothing the 
spectra. 

Annihilation modes involving quarks, heavy vectors or the Higgs are disfavored. A precise 
observation of the excess around its end-point at ~ 2 TeV would allow to settle this issue 
and to partially disentangle the allowed leptonic modes. Additionally for all models that are 
consistent, the PAMELA positron fraction should continue to increase over the entire range of 
energy that PAMELA can explore. 

We find that, in view of the FERMI observation, the energy spectrum of photons due to 
Inverse Compton scattering of can now be robustly predicted. FERMI should be able to see 
a 7 excess if the are generated by DM in the DM halo, rather than locally (e.g. by a nearby 
pulsar). This allows one in principle to probe the DM interpretation of the excesses and 
possibly learn about the DM density profile. To extract further information about the particle 
physics properties of DM that can explain the excesses one would need to observe or more 
strongly bound the associated neutrino and final state radiation 7 fiuxes. 
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Figure 13: Function f^{Ei,) relevant for the neutrino signals. 



A Neutrinos from the Galactic Center 

SK observes the flux of through going (TG) muons, generated by z/^ scatterings in the matter 
below the SK detector down to a depth given by the muon range. We recall that muon energy 
losses and range R in matter can be approximated as 

^ = a, + /3,E„ R{E,) = Un{l + (^E,) (15) 
ax p 

where = 2 ■ 10^^ GeVg^^ cm^ and = 3 ■ lO^^g^^ cm^ becomes dominant at E^^ > TeV. 
The h'fj,N — )■ fiN' cross section is 



da 2Gpmjv 



dE^ 71 



Pi +P2 



E,, 



(16) 



where pi^2 describe the parton content of the nucleon: pi 0.2 and p2 ~ 0.05 for z/, and 
the reverse for i/, having assumed an average nucleon N with equal number of n and p and 
Ey <^ M^/rriN ~ 10 TeV. Then, the total number of TG muons produced by a generic energy 
distribution dNy/dEy of neutrinos is: 



= / dE^^ii^fiEy) (17) 



dN, 







dE, 



V 



where the dimensionless function, /, is given by 

''^^ . „ da 



f{E,) = A / " dE^-^[R{E^) - i?(Etre.h)] 

and can be written as, 



m) = E f'-^f'+P^h n{E.). (19) 

Here //j = 1 as long as i?^ ^ TeV such that one can approximate = 0. We plot the function 
//J in Fig. 13. It vanishes at negligibly small neutrino energies due to the (negligibly small) 
threshold -Etresh = 1-6 GeV in SK. 
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We therefore find, 




TG 



2 47rM2 



JSf X 



M2G|(2pi+p2) 



27ra 



Sf 



1 



dx X 



^dN, 



dx 



V 



(20) 



Notice that the dependence on M cancels out, up to the residual dependence implicit in the 
factor //J. Up to this factor our formulae confirm the ones in [38]. 

B Details of the Fitting Procedure 

We smoothly vary between the NFW, Einasto {a = 0.17) and isothermal DM density profiles, 
defined as in [11]. In particular, we keep fixed the local DM density Pq = 0.3 GeV/ cm^, which 
is uncertain by about a factor of 2. Since the indirect signals are proportional to (o"f)p^ or 
to p/t for DM annihilations and decays respectively, different values of can be studied by 
rescaling {av) or r. 

We approximate the diffusion zone with a cylinder of half-thickness L centered on the 
galactic plane and adjust the e^,p diffusion parameters as function of L in order to reproduce 
cosmic ray data. We smoothly vary between the favored case L = 4 kpc (MED configuration) 
and the extremal cases L = 1 kpc (MIN configuration) and L = 15 kpc (MAX configuration). 
All details are as in [2]. 

Concerning the astrophysical backgrounds, we fit in two ways. First in the method similar 
to [2], which employs the background estimate from [53], allowing for an uncertainty in the 
overall normalization and the spectral index of each component. The new FERMI data allow 
us to now enlarge the uncertainties, as the fit now relies more on data rather than the theoretical 
computations. Second we use an almost completely data driven approach to backgrounds as 
described in [12]. We find that for all fits these two methods give similar results and we consider 
our background estimates to be as robust as possible without the ability to determine a better 
estimate of backgrounds independently. 

We checked that including the latest 2009 HESS e'^ + e' [4] data negligibly affects our fits, as 
the latest HESS data extend to lower energies where FERMI provides more precise observations. 
Similarly, the preliminary PAMELA measurement of the un-normalized e~ flux [17] similarly 
has a negligible impact on the global fits, in view of the more precise FERMI measurements. 
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